Dynamic ordering of driven vortex matter in the peak effect regime of amorphous 

MoGe films and 2H-NbSe 2 crystals 
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Dynamic ordering of driven vortex matter has been investigated in the peak effect regime of 
both amorphous MoGe films and 2H-NbSe2 crystals by mode locking (ML) and dc transport mea- 
surements. ML features allow us to trace how the shear rigidity of driven vortices evolves with the 
average velocity. Determining the onset of ML resonance in different magnetic fields/temperatures, 
we find that the dynamic ordering frequency (velocity) exhibits a striking divergence in the higher 
part of the peak effect regime. Interestingly, this phenomenon is accompanied by a pronounced 
peak of dynamic critical current. Mapping out field-temperature phase diagrams, we find that di- 
vergent points follow well the thermodynamic melting curve of the ideal vortex lattice over wide 
field/temperature ranges. These findings provide a link between the dynamic and static melting 
phenomena which can be distinguished from the disorder induced peak effect. 



PACS numbers: 74.25.Qt, 74.25.Sv, 74.25.Fy 



1. INTRODUCTION 



Vortex matter in type II superconductors has been rec- 
ognized as an ideal system for studying collective trans- 
ports of periodic manifolds driven over the pinning en- 
vironments [1] [2]. In particular, much recent attention 
has focused on the issue of dynamic ordering, separat- 
ing elasticity dominated, coherent lattice flow at large 
velocity, from plastically deformed, incoherent defective 
flow at small velocity. Theoretical considerations for this 
phenomenon were first proposed by Koshelev and Vi- 
nokur by introducing a concept of shaking temperature, 
which characterizes fluctuating motion of vortices due to 
the interaction with disorder quenched in a host mate- 
rial [3]. Assuming the shaking temperature to decay as 
ksTsh = r/v on increasing the average velocity v, they 
proposed that the ordering (crystallization) velocity di- 
verges as v c oc r/fc_s(TV — T) with a disorder parameter 
r on approaching the thermal melting temperature Tm 
from below. Later, more sophisticated theoretical [4] [5] 
[6] [7] and numerical studies [8] [9] [10] [11] [12] [13] have 
suggested various coherent flow states with different po- 
sitional and orientational order like moving glass state 
and smectic flow state, and some of them have proposed 
multiple transitions in ordering process. 

While experimental evidences for the ordering were re- 
ported in neutron diffraction and Bitter decoration stud- 
ies [14] [15] [16] [17], investigations of the divergence 
at the ordering point in dc transport studies were only 
based on the assumption that a S-shaped anomaly in 
dc current- voltage characteristics marks current-induced 
dynamic ordering of driven vortices. In NbSe2 crystals 
[1] [2], the inflection anomaly is observed in the peak 
effect regime, where the critical current shows a peak in 
the vicinity of the second critical field H C 2- Meanwhile in 
amorphous (a-)MoGe thin films [18], instead of the peak 



effect, the anomaly is accompanied by a peak in the "dy- 
namic" critical current, which characterizes a change of 
dynamic friction of vortices in the flux-flow state. These 
studies claimed the divergence of the ordering current 
around these peak behaviors and suggested the relevance 
of the thermal melting to the peaks of the static or dy- 
namic pinning force. However, a study of scanning Hall 
microscopy has shown that geometrical inhomogeneity 
in sample edges leads to a macroscopic coexistence of 
ordered and disordered flow states and their dynamic 
evolutions with transport current may result in similar 
S-shaped current- voltage characteristics [19] [20]. Thus, 
it has remained controversial if the S anomaly marks the 
dynamic order of driven vortices. 

Recently, mode locking (ML) has proved to be a pow- 
erful experimental technique for studying the dynamic 
ordering of driven vortex matter. ML is a dynamic res- 
onance between internal lattice modes of a driven lat- 
tice and an ac drive superimposed on top of dc drive 
[21] [22] [23] [24] [25] [26] [27]. This technique directly 
probes the shear rigidity of driven vortices as a function 
of the velocity and allows us to find the dynamic ordering 
(or melting) point from the onset (or disappearance) of 
ML resonance. As demonstrated in a recent ML study 
on NbSe2 crystals, the ordering velocity exhibits a steep 
increase in the peak effect regime and this defines un- 
ambiguously the thermal melting point for driven vortex 
lattice [27]. 

To elucidate the relevance of thermal melting to the 
peak behaviors of the critical current and the dynamic 
critical current, we will present in this article a system- 
atic study of ML and dc transport measurements on both 
cv-MoGe films and NbSe2 crystals over wide magnetic 
field/temperature ranges. We will show how the ordering 
frequency (velocity) grows with field and diverges in the 
peak effect regime. We will map out field-temperature 



1 



phase diagrams and show how the divergent field, to- 
gether with the peak fields of the critical current and 
the dynamic critical current, grow with lowering temper- 
ature. We will discuss influences of quenched disorder 
and thermal fluctuations on the onset of ML resonance 
and compare the divergent points to the thermodynamic 
melting curve of the ideal vortex lattice proposed by a 
quantitative Ginzburg-Landau (GL) model [28]. A brief 
comment on the peak behavior of the dynamic critical 
current will be made. 



2. EXPERIMENTAL 

Amorphous Moi-aGe^ films (x — 0.22 ± 0.01%) were 
sputtered on silicon substrates mounted on a rotating 
copper stage held at room temperature by water cool- 
ing. The films were structured into a Hall-bar shape by 
a liftoff technique. After partly deposition of silver on top 
of the films, thin gold wires were connected to the films 
by using small indium pieces. We used two films (MGjjl 
and MG|j2) with nearly identical geometry (1.2 mm in 
length I, 0.3 mm in width w, and 0.33 /tm in thickness 
d) and material parameters of superconducting transi- 
tion temperature T c , normal resistivity p n and a slope 
of the second critical field near T c S(= p^dH^/dT \t c ), 
which are given in table 1 . We determined T c from resis- 
tive transition by Aslamazov-Larkin fluctuation theory 
[29]. H c2 for magnetic field perpendicular to the films 
was determined from the lowest-landau-level (LLL) scal- 
ing analysis of fluctuation conductivity [30] [31]. The 
average composition for the films was obtained by the 
electron probe micro analyzer. 

We also used pure 2H-NbSe 2 single crystals grown by 
an iodine vapor transport method [32]. The thinnest 
platelets with optically flat surfaces on both sides were 
cut in a Bar-shape. We glued these on silicon substrates 
with negative resist and cleaved them without observable 
surface damage. After deposition of silver partly on crys- 
tal surface, thin gold wires were connected to voltage and 
current contacts by silver past and indium solder, respec- 
tively. The geometry and material parameters for two 
cleaved crystals (NSjJl and NSJJ2) used in this study were 
also summarized in table 1. T c was determined by the 
mid-point of the resistive transition, which has a transi- 
tion width of about 50 mK between 10% and 90% of the 
normal state resistance at 8 K, just above T c . Residual 
resistance ratio defined by i?(295 K)/i?(8 K) is about 30. 
We applied magnetic field perpendicular to the ab plane 
of the crystals. H C 2 was determined from the intersection 
between the linear extrapolation of flow resistance and 
the normal resistance [33] . The thickness of the cleaved 
crystals was estimated from the room temperature resis- 
tance by d — (7/u>)/9 n (295K)/i?(295K) using the room 
temperature resistivity p n (295K) = 1.2 x 10~ 6 £!m deter- 
mined before [27]. 

We used Oxford cryostat with a VTI inset and a 17/19 



superconducting magnet. We glued a sample on a cop- 
per stage at the bottom of our measurement insert, and 
measured temperature by using a calibrated Cernox sen- 
sor attached to the stage. For measurements below 4.2 K 
the sample was immersed in liquid of 4 He introduced in 
the VTI insert. We regulated the vapor pressure of the 
liquid by using a throttle valve, a capacitance manome- 
ter, a PID controller and a rotary pump. Above 4.2 K, 
we used a temperature controller (Cryocon 62), which 
regulates temperature of the vaporizer at the bottom of 
the VTI insert and that of the sample stage of the mea- 
surement insert. We employed four-terminal pair method 
with a precision LCR meter (HP 4285A/001) and nano- 
voltmeters (Keithley 2182). 

The sample was wired with four-coaxial cables of which 
guards are connected. For a precise measurement of the 
ac current in the frequency range of 0.1-30 MHz, we cal- 
ibrated the coaxial cables in two steps: First we made 
open and load calibrations by unwiring and wiring a 100 
ohm chip resistor to the cables at the sample stage, re- 
spectively. Then, the sample was connected to the cables 
and cooled to the lowest temperature of ~ 1.2 K, which 
is far below T c of a-MoGe films and NbSe2 crystals used 
in this study. Ensuring negligible contact resistance and 
undetectable resistance of the sample, we made short cal- 
ibration as function of frequency as the second step. 

During the ML measurement, we probed the dc volt- 
age by ramping the dc-current up or down with a super- 
imposed ac current on top. The dc current value was 
obtained by measuring the dc voltage over a standard 
resistor. Other dc transport measurements of resistive 
transition and dc current-voltage characteristics without 
ac current were performed by using a dc current source 
(Keithley 220) and nanovoltmeters with proper filtering. 



3. RESULTS AND DISCUSSIONS 

Figure 1 displays a plot of the critical current I c vs 
field measured in 1.96 K for MGjjl. Here I c is deter- 
mined from dc current-voltage (I — V) characteristics by 
a voltage criterion of 1 pV. As observed, on increasing 
magnetic field I c starts to increase around 5.5 T (the 
onset field) and exhibits a broad peak around 7 T (= 
jjL{)H p ). Above 7 T, it drops rapidly and vanishes around 
7.8 T. We also plot the dynamic critical current I c ,dyn, 
obtained by extrapolating the flux-flow behavior linearly 
to the zero voltage (for definition, see the inset to Fig. 
1). This characterizes the pinning strength (or dynamic 
friction) for driven vortices in the flux-flow state. As ob- 
served, I c ,dyn exhibits a pronounced peak against field: 
It starts to show a rapid increase around 7 T and ob- 
tains a sharp maximum at about 7.6 T (= p.oH Pt d yn ). 
Interestingly, the peak field (onset field) of the dynamic 
critical current does not coincide with the corresponding 
peak (onset) field for the critical current. The former ap- 
pears in the higher part of the peak effect regime of the 
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latter, i.e., H p .dyn > Hp- Such peak anomaly of the dy- 
namic critical current is also observed in NbSe2 crystals. 
In the conventional picture, the peak in the critical cur- 
rent marks the pinning induced, structural transforma- 
tion from weakly disordered to amorphous arrangement 
of vortices [45] [36] [37]. Thus, from a phcnomcnologi- 
cal point of view, another structural transformation in 
the flux-flow state may occur at the peak in the dynamic 
critical current. 



Mode Locking 

In order to explore the dynamic structural change in 
the flux-flow state, we have employed the ML technique. 
To show typical ML resonant features, we first present 
ML results below the peak effect regime. Shown in Fig. 
2 is a series of I — V curves of MGj}2 measured with 
superimposed 10 MHz ac current of various amplitudes 
at 4.2 K and a field of 2.0 T (< ^ # P (4.2K) = 3.0 T). 
A curve denoted as DC represents the pure dc I — V 
curve without ac current I ac . The application of I ac in- 
duces clear ML steps at equidistant voltages denoted by 
p/q=\/l and 2/1. The first step corresponds to the fun- 
damental (p/q =1/1), and the other is a higher harmonic 
(p/q=2/l). Also even higher harmonics are observed at 
higher voltages, but are not shown. Because of the ML 
resonant condition of v = (p/q) fa with integers p, q and 
the lattice spacing of vortex lattice a, these ML volt- 
ages are directly proportional to the resonant frequency 
of ac current according to V p / q = vBl = (p/q)$ fl/a± 
with the row spacing of vortex lattice a±(= <fr /aB) and 
the vortex density B [27]. Step width at the funda- 
mental (and also harmonic) ML resonance varies with 
I ac , exhibiting a oscillatory behavior like squared Bessel 
function of the first kind of the first (higher) order [26] . 
Since the oscillatory period depends on frequency, tem- 
perature, and magnetic field, we use for consistency in 
the following analysis the maximum width AI max of the 
fundamental resonance as the amplitude of the ML step 
width. 

In Fig. 3(a) we show how AI max varies with frequency. 
At lower frequencies (< 10 MHz), the ML step width in- 
creases linearly with /, while at higher frequencies (> 10 
MHz) it levels off and becomes independent of /. Such 
behavior is reasonably described by an empirical func- 
tion of AI m ax = / s tanh(//0.7/ p ), displayed by the solid 
curve [26]. From the fit we obtain a pinning frequency 
f p = 13 MHz and a saturated current width I s — 16/iA 
[38]. However, the ML width at lower frequencies de- 
viates downward from the empirical function. A linear 
extrapolation of the data indicated by the dotted line 
seems to reveal a threshold behavior of the ML width, 
i.e., there is a finite frequency (velocity) above which the 
ML resonance can be observed. 

Such threshold behavior is pronounced in the peak ef- 
fect regime, as can be seen in Fig. 3 (b), where the fre- 



quency dependence is shown of the ML width measured 
at 4.2 K and magnetic fields of 3.25 T, 3.4 T, and 3.45 T 
(above [ioH p = 3.0 T). For instance, the linear extrapo- 
lation at 3.25 T clearly displays the threshold behavior at 
an onset frequency f c w 1.4 MHz denoted by an arrow. 
Below f c , no ML resonance feature is detectable at any of 
the ac current amplitudes, indicating the absence of shear 
rigidity in the driven vortex matter. Thus, f c marks a 
dynamic ordering frequency (or velocity) [25] [26]. When 
field is increased by 0.15 T to 3.4 T, f c increases by 1.7 
MHz to w 3.1 MHz. With further increase of field by 
only 0.05 T to 3.45 T, f c shows a rapid increase by s=s 2.2 
MHz and becomes ~ 5.3 MHz. In slightly higher fields 
an even more rapid increase of f c is observed. 

We summarize results of the onset frequency at 4.2 K 
in a semi-logarithmic plot of f c vs. magnetic field, dis- 
played in Fig. 4(b). In low fields, f c (depicted by solid 
diamonds) increases slowly with field. Around the peak 
field of the critical current (=3.0 T), it turns to a rapid 
increase, followed by a small knee around 3.2 T. After 
the knee, it exhibits a diverging behavior towards just 
above 3.5 T. This is also displayed in a linear plot given 
in the inset of Fig. 4(b). The open symbols represent 
the onset frequency determined from the magnetic field 
dependence of the ML width [27]. Both data sets show 
the same diverging behavior of f c (H). 

The divergence of f c is in qualitative agreement with 
the dynamic ordering picture proposed by Koshelev and 
Vinokur [3]. The dynamic crystallization frequency can 
be given by [25] [26] 
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with a disorder parameter j u and the flux-flow resistiv- 
ity pf. This expression qualitatively explains the diver- 
gence of f c (or v c ) against temperature observed in an 
a-MoGe film [39] and also in mesoscopic flow channels 
[25]. Moreover, as displayed by the solid curve in Fig. 
4(b) (also in the inset), we find that the diverging be- 
havior of f c against magnetic field is approximated well 
by a similar diverging function of f c = fo/(l — H/Hm) 
with HoHm — 3.53 T (indicated by the broken vertical 
line) and fo = 0.11 MHz over more than one decade in 
frequency range. These qualitative agreements allow us 
to identify the field (or temperature) of the divergence 
as the thermal melting point of the coherently flowing 
vortex lattice [27]. Away from Hm at lower fields, the 
situation is different. Here the onset of ML resonance 
occurs at small frequencies (velocities) and the influence 
of quenched disorder dominates the ordering transition. 
The contribution of quenched disorder on f c may be de- 
termined from following arguments: The disorder param- 
eter is given by the amplitude of the pinning correlation 
function of the random pinning potential which decays 
over the pinning radius r p [3] . Using the maximum pin- 
ning force F p per pinning site, we may write 
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where N p (= a 2 dn p ) is the number of pinning sites per 
vortex line volume, n p is the density of pinning cen- 
ters, W(= ripFp/2) is the Larkin-Ovchinnikov pinning 
parameter [40] and r p a/2 [41]. Assuming Tm — T w 
Po(H M — H)/S, we find 



fc ~ /o 



with 



/o = 



1 



(1 - H/H M ) 
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(3) 



(4) 
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We know that the pinning properties for a-MoGe films 
are well described by the delta- T c pinning, i.e., W — 
Cb(l - b) 2 with a reduced field b(= H/H c2 {T)) [41] 
[44] [45]. The flux flow resistivity for the amorphous 
films in the GL regime can be written as {pf/pn)^ 1 = 
F{b)/\Jl — t + 1 with a reduced temperature t = T/T c 
[34]. In fields measured here, i.e., b = p n H / ' ^H c2 > 1.4 
T/4.54 T, the function F(b) can be approximated by an 
interpolation formula F(b) = (1 - fe) 3 / 2 (0.43 + 0.69(1 - 
&))/6[43]. Substituting these formula in Eq. (4), together 
with poH = B, we find that f c varies with magnetic field 
as 

(1-b) 2 1 



fc OC 



(5) 
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As displayed by the dotted curve, this shows a slow in- 
crease with field away from H ss Hm and explains our 
results below H p . The proportionality constant extracted 
from the fit provides a rough estimate for the coefficient 
C in the expression for the pinning strength W , namely 
C w 1.2*10- 8 N 2 /ni 3 at T w 0.7 T c . This value is reason- 
ably consistent with values obtained for previous trans- 
port measurements on similar amorphous superconduct- 
ing films [41] [44] [45]. Just above H p , f c switches to a 
rapid increase up to the knee. This behavior follows nei- 
ther the diverging behavior close to Hm displayed by a 
solid curve nor the low field behavior by a dotted curve. 
It indicates a rapid enhancement of the disorder contri- 
bution (probably the pinning parameter W). Since the 
ordering frequency (velocity) around H p is small and the 
disorder effect does not diminish significantly, this jump 
behavior could be related to the underlying mechanism 
for the peak effect, i.e., the enhancement of the pinning 
force caused by structural disordering transformation of 
static vortices at H p [45] [36] [37]. 

We have systematically made ML and dc transport 
measurements in different temperatures on both a-MoGe 
films and NbSc2 crystals. Results commonly show that f c 
exhibits the jump and subsequent knee behavior around 
H p . The diverging behavior of f c occurs always in the 
higher part of the peak effect regime and the divergent 
point lies in between the peak field and the second criti- 
cal field, i.e. H p < Hm < H c2 . Thus, observed behavior 
of f c allows us to separate the thermal melting marked 
by the divergent point from the disorder induced trans- 
formation marked by H p . 



Phase diagrams 

Based on these findings, let us map out a field- 
temperature phase diagram for a-MoGe films of MGjjl 
and MGjj2. To show results in two films together in Fig. 
5(a), magnetic field and temperature are normalized by 
the GL second critical field Hg L (0)(= ST C ) and T c , re- 
spectively. For clarity, results of MGjjl and MGj)2 are 
depicted by black and red colored symbols, respectively. 
Good agreement is visible in temperature dependencies 
of h c2 (= H c2 /Hf 2 L {Q)) plotted for the films with solid red 
and black circles: On lowering temperature from t = 1 
h c2 (t) increases linearly with 1 — t displayed by a bro- 
ken line. For t < 0.8, h c2 (t) deviates downward from 
the linear dependence. Such behavior is well described 
by the mean field curve (displayed by the solid curve) for 
the dirty limit in Werthamer-Hclfand-Hohenbcrg (WHH) 
theory [48], which is given by 



ln(t) = *(l/2) - *[l/2+ (2/ir 2 )h c2 (t)/t] 



(6) 



where $ is the digamma function. 

In Fig. 5(a) the results of h M (t)(= H M /Hg L (0)) for 
the films are also depicted by solid squares. Focusing 
on results of MGjjl, we find that on lowering tempera- 
ture from t = 1 Km (t) increases slowly and field distance 
between /im(0 and h c2 (t) becomes larger. With further 
lowering t, fiM(t) gradually increases faster and it be- 
comes nearly parallel to h c2 (t) for t < 0.8. Thus, the 
field distance between /im(*) and h c2 (t) is nearly con- 
stant for t < 0.8. Results of MGjJ2 follow the h M (t) 
curve for MGjjl well. Because of the thermal melting of 
the flowing vortex lattice marked by liM{t), the regime 
in the phase diagram between /im(^) and h c2 (t) should 
represent a flowing liquid of driven vortices. 

We also map out the ficld-tempcraturc phase diagram 
for NbSc 2 crystals of NSjjl and NSjj2 in Fig. 5(b). h c2 (t) 
increases linearly with 1 — t down to the lowest temper- 
ature measured. This is consistent with thermodynamic 
studies on NbSe 2 crystals [46] [47]. In contrast to the 
results on a-MoGe films, tiM(t) lies just below the mean 
field line except for low temperatures (t < 0.4). For in- 
stance, the field distance between /im(^) and h c2 (t) ait^ 
0.7 is 0.02 for NbSe2 crystals, which is smaller than 0.08 
for a-MoGe films. In other words, the vortex liquid flow 
regime for NbSe2 crystals is much narrower than that for 
a-MoGe films. 

This difference should be related to the strength of 
thermal fluctuations, which may be characterized by the 
Ginzburg number given by 



Gi = 



1 



2 V4^o^c(0) 2 ^ h (0) 



(7) 



where 7 is an anisotropy parameter, H c (0) is the thermo- 
dynamic critical field at T = 0, and £ a b(0) is the coher- 
ence length for the ab plane at T — [42] . For estimation 
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of Gi on a-MoGe films, we use the dirty-limit expression 
for mo# c (0) = ST c /2.54k with k = 3.54 * l0 4 [p n S} 1/2 
[41] and 7=1. For NbSe2 crystals, we use values of the 
GL thermodynamic field noH c (0) = 0.23 T and 7 = 3, 
reported in the literature [46] [47] . As listed in Table 1 , a- 
MoGc films have a nearly identical value of Gi ss 8 * 10~ 5 
which is two orders of magnitude larger than that for 
NbSe2 crystals. This is qualitatively consistent with the 
broader regime of the vortex liquid flow for a-MoGc films. 

Quantitative estimation of the distance between the 
thermodynamic melting curve and the mean field line 
can be made by using the recent quantitative GL result 
proposed by Li and Rosenstein (LR) [28] . In their theory, 
the thermodynamic melting transition sets in when the 
LLL-scaled temperature given by 



«t = -a h yy -jp-nthj (8) 

is equal to -9.5. Here Ng%(= Gi/ir 2 ) is their Ginzburg 
number and a^(= 1 — t — h) is the distance from the mean 
field line. Substituting Gi for a-MoGe films and NbSe 2 
crystals from table 1, we calculated how far the melting 
curve is located from the mean field line (1 — t — h C 2(t) = 
0). A solid blue curve in Fig. 5(a) represents the cal- 
culated melting curve for MGjJl. Good agreement with 
h,M(t) data is visible at high temperatures (t > 0.6). We 
note that the melting curve at low temperatures (t < 0.6) 
is not shown since clear deviations of h C 2 from the linear 
dependence, which is presumed by the model, appear at 
low temperatures. We also display the calculated melt- 
ing curve for NSjjl in Fig. 5(b). Because of small Gi, the 
calculated melting curve lies just below the h C 2(t) line. 
Again, good agreement of Km (t) with the melting curve 
is visible at high temperatures (t > 0.6). Taking into 
account that this model is valid in temperatures not far 
from t = 1, the agreements of /im(£) with the LR model 
for both MGjtl and NSjJl having more than two orders 
of magnitude difference in Gi are quite remarkable, pro- 
viding a further link between the thermal melting point 
of the flowing lattice and the thermodynamic one of the 
static and ideal vortex lattice. We note that the influence 
of applied currents upon the thermodynamic properties 
is negligibly small [49]. 

We also plot the peak fields h p (= H p /H^ 2 L (0)) of the 
critical current (open circles) in Fig. 5(a) and (b). In 
both cases they lie below the thermodynamic melting 
curve (or IiM(t))- This indicates that for static vor- 
tices the pinning induced, structural transformation into 
amorphous array of vortices occurs prior to the thermo- 
dynamic melting and therefore a crossover (or a contin- 
uous transition) from glassy to liquid states may occur 
around the thermodynamic melting point. We empha- 
size that in the higher part of the peak effect regime the 
dynamic melting features discussed above is observable 
only when vortices are driven fast enough and that the 
elasticity is recovered by the dynamical ordering from the 
amorphous static state. 



Finally we comment on the peak behavior in the dy- 
namic critical current. We find that all the peak points of 
the dynamic critical current observed for both a-MoGe 
films and NbSe2 crystals coincide nearly with correspond- 
ing divergent points of ftjw(i) ( n °t shown for clarity). 
This agreement is exemplified in Figs. 4(a) and (b), 
where the peak point of the dynamic critical current oc- 
curs very close to the divergence of f c . These coinci- 
dences suggest that the peak point of the dynamic crit- 
ical current could be a good indication for the thermal 
melting of the coherent lattice flow. However, it is not 
clear to us why the lower part of the peak anomaly ob- 
served in both a-MoGe films and NbSe2 crystals change 
gradually and do not show a jump at the thermal (or 
thermodynamic) melting point. This gradual change is 
possibly be related to an extrinsic origin like inhomogene- 
ity in samples and/or geometrical influences of edge [51] 
and surface roughness [52]. 



SUMMARY 

In summary, employing the mode locking technique, 
we have presented the dynamic ordering of driven vor- 
tex matter in the peak effect regime of both amorphous 
MoGe films and NbSe2 crystals. The dynamic ordering 
frequency (velocity) marked by the onset of the ML res- 
onance exhibits a sharp increase and subsequent knee 
behavior around the peak field of the critical current, 
followed by a striking diverging behavior in the higher 
part of the peak effect regime. The former jump behav- 
ior indicates the rapid enhancement of pinning contri- 
bution, probably originating from the underlying struc- 
tural transformation of the static vortex lattice due to 
quenched disorder at H p . Meanwhile the latter diverging 
behavior is in good agreement with the dynamic order- 
ing model and the divergence point is identified as the 
thermal melting point of the driven vortex lattice. 

The thermal melting point measured at each tem- 
perature lies in between the peak field and the sec- 
ond critical field, and it coincides nearly with the peak 
point of the dynamic critical current. Mapping out the 
field-temperature phase diagrams for amorphous MoGe 
films and NbSe2 crystals, we find that the thermal 
melting points follow well the thermodynamic melting 
curve for the static and ideal vortex lattice over wide 
field/temperature ranges. These findings clearly sepa- 
rate the peak effect from thermal melting and provide a 
link between static and dynamic melting phenomena. 
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FIGURE CAPTION 



Fig. 1 Magnetic field dependence at 1.96 K of criti- 
cal current I c (solid circles) and dynamic critical current 
I c ,dyn (open circles) for an a-MoGe film (MGjjl). In the 
inset an / — V curve measured at 7.2 T is given. The 
solid line represents the linear extrapolation of the flux 
flow behavior to the zero voltage. As indicated by an 
arrow, the intersectional point defines I c .dyn- 

Fig. 2 A set of dc current-voltage curves for an a- 
MoGe film (MGjj2) measured at 2.0 T in 4.2 K with 
superimposed ac currents (amplitude 3.3, 2.6, 2.1, 1.6, 
1.3, 1.0, 0.82, 0.65, 0.52, 0.41, 0.33, and 0.26 mA from 
left to right). A curve denoted as DC indicates the pure 
dc I — V curve. Application of ac current induces clear 
mode-locking current steps at equidistant voltages de- 
noted by p/q =1/1 and 1/2. 

Fig. 3 Frequency dependence of maximum current 
width AI max of the fundamental ML step for an a-MoGc 
film (MG|J2) measured (a) at 2.0 T, far below the peak 
effect regime, and (b) at fields of 3.25 T, 3.40 T and 3.45 
T in the peak effect regime in 4.2 K. A solid curve rep- 
resents an empirical function for the ML width given in 
text [26] . Dotted lines represent the linear low frequency 
behavior of the ML width. Arrows mark the onset fre- 
quency f c for ML resonance. 

Fig. 4 Magnetic field dependence of(a) critical current 
I c , dynamic critical current I c ,dy n and (b) the onset fre- 
quency f c for an a-MoGc film (MGjj2) at 4.2 K. Peak 
points of the critical current and the dynamic critical 
current are indicated by arrows in (a). The inset to (b) 
displays the plot of f c on a linear frequency scale. 

Fig. 5. (Color online) Magnetic field-temperature 
phase diagrams for (a) a-MoGe films and (b) NbSe2 crys- 
tals. The magnetic field and temperature are normalized 
by superconducting transition temperature T c and the 
Ginzburg-Landau second critical field H^ L (0)(= ST C ). 
Plotted are reduced second critical fields h c2 (small solid 
circles), divergence fields hu of the ordering frequency 
(solid squares) and peak fields h p of the critical current 
(open circles). In (a), black and red colored symbols cor- 
respond to results for MGjjl and MGjj2, respectively. In 
(b), black (red) colored symbols are for NSjjl (NSjj2). In 
(a) a solid black curve represents the mean field curve for 
the dirty- limit in WHH theory [48]. For comparison, the 
linear dependence of h c2 on 1 — t is displayed by a broken 
line in (a) (also in (b)). A dotted curve in (a) is guide to 
the eye. Blue curves in (a) and (b) represent thermody- 
namic melting curves obtained from a quantitative GL 
model [28] for MGjjl and NSjjl crystals, respectively. 

Table 1. Properties of amorphous MoGe films and 
NbSe2 crystals. 
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TABLE I. Properties of amorphous MoGe films and NbSe2 
crystals 



Material Sample 


T C (K) Pn - 


(8 K)(fj,Qm) 


-HodH c2 /dT\ Tc (T/K) !(mm) 


ra(mm 


) d(nm) 


Gi 


MoGe MGfll 


6.11 


1.84 


2.6 


1.2 


0.3 


0.33 


8.9*10 


MGP 


6.05 


1.72 


2.7 


1.2 


0.3 


0.33 


8.1*10 


NbSe 2 NSttl 


7.1 


0.040 


0.75 


0.8 


0.4 


1.8 a 


6.1*10 


NSJJ2 


7.2 


0.034 


0.72 


0.72 


0.63 


0.9 a 


5.8*10 



a The thickness was estimated from the normal resis- 
tance at the room temperature _R(295K) by d — 
(Z/w)p„(295K)/.R(295K) with p„(295K) = 1.2/ifim [27]. 
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